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Abstract 

We present a measurement of the CP violation parameters in Kgir^ir^ and —>■ KgKg 

decays using a data sample containing 657 x 10^ BB pairs collected with the Belle detector at the 
KEKB asymmetric-energy e+e~ collider operating at the T{4S) resonance. We measure 

^ROnOnO = +0.43 ± 0.49 ± 0.09, 

^^o^o^o = -0.17 ±0.24 ±0.06, 

Skoko = -0.38 ± 0.77 ± 0.08 and, 

AkOk° = -0.38 ± 0.38 ± 0.05, 

where the first and second errors are statistical and systematic, respectively. 
PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



3 



In the standard model (SM) framework, CP violation arises only from the irreducible 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix f]]]. CP violation in flavor chang- 
ing neutral currents (FCNC) is sensitive to new physics (NP), as new particles beyond the 
SM may be involved in loop diagrams. Such additional loop diagrams potentially add new 
CP phases and induce deviations from the SM expectation for time-dependent CP asym- 
metries jH]. 

In the decay chain T(4S') B^B^, one of the two B mesons decays into a CP eigenstate 
fcp at time tcp, and the other decays into a flavor specific state /tag at time ttag- The CP 
violation parameters are measured using the time-dependent decay rate ^ 

g-|At|/rgO 

P(At) = [1 + g • {5/ sin( AmrfAt) 

+AfCos{AmdAt)}], (1) 

where At = tcp — ttag, t^o is the B^ lifetime, Arrid is the B^B^ mixing frequency and q = +1 
(— 1) when /tag = -B°(i?°). The parameters Sf and Af represent mixing-induced and direct 
CP- violation, respectively. 

In the SM, the CP violation parameters in 6 — s "penguin" and b ^ c "tree" transitions 
are predicted to be Sf ~ — ^/sin20i and Af ^ with small theoretical uncertainties 
where 0i is one of the CKM weak angles and ^/ is the CP eigenvalue of the final state. 
Recent measurements [sl, 0, 0| however, indicate that the effective sin20i value, sin2(j)'f, 
measured with penguin processes is different from sin20i = 0.687 ± 0.025 measured in tree 
decays by 2.6 standard deviations (cr) [8i]. New particles in the loop diagrams may have 
shifted the weak phase. 

Belle and BaBar recently observed b ^ d FCNC transitions in B meson decays into 
two kaons B+ K+K^ and B^ K^K^ i, 0, El. The CP violation parameters in 



b ^ d transitions are expected to be small because the weak phases in B^B'^ mixing and 
the b ^ d transition cancel [l2|. Thus a measurement of large CP violation could indicate 
the existence of NP. It has been also pointed out that the CP violation parameters as well 



as the branching fraction for B^ — > K^K^ decay can be used to constrain 02 131 . 

In this report, we present measurements oi CP violation parameters in S° — >■ K^t^^tt^ 
decays that proceed through b — > sqq [q = u, d) transitions, and B^ — *• K^K'g decays in 
which pure b dss penguin processes are the dominant contribution. Because the ^/ for 
B^ — > K^ir^ir^ decays is +1 for any intermediate resonance IJ], the SM expectations for the 



CP violation parameters are ^ ~ and S = — sin 20i with small theoretical uncertainty |15| . 
Both analyses employ a technique [l6|] that reconstructs the B^ decay vertex position using 
only Kg^s decaying into 7r"'"7r~. The CP violation parameters of these B decay modes have 
been measured by BaBar 0, [ij with 227 (348) x 10*^ BB pairs for S° ^ i^gTr^vrO (i^^K^). 



They found a 2.2a deviation from the SM expectation for B^ K^tt^tt^ and a large iS^o^^o 
value; however, no clear conclusion about NP can be obtained because of the large statistical 



errors. 



We use a data sample containing 657 x lO^BB pairs collected by the Belle detector 
at the KEKB e^e~ asymmetric-energy (3.5 on 8 GeV) collider operating at T(4S') reso- 
nance ISj. The T(45') is produced with a Lorentz boost factor of j3'~f = 0.425 approximately 
along the electron beamline {z). 

Since the two B mesons are produced nearly at rest in the T(45') center-of-mass system 
(cms), we determine At = Az/jSjc, where Az is the distance between the two B decay 
positions in the z direction. In the Belle detector, a silicon vertex detector (SVD) and a 
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50-layer central drift chamber (CDC) are used for charged particle tracking. Photons are 
detected with an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals. The 
devices are placed inside a superconducting solenoid coil providing a 1.5 T magnetic field. 

We reconstruct — > tt+vt^ candidates from pairs of oppositely charged tracks. The 
candidates are required to satisfy the following conditions: (a) the two charged tracks do 
not originate from the beam interaction point (IP), (b) the z positions of the two charged 
tracks at the K^g decay position have to match, (c) the direction of the pion pair momentum 
is consistent with the direction from the IP to the decay position, and (d) the invariant 
mass of the two charged tracks is within ±15 MeV/c^ of the nominal mass. The criteria 
(a) to (c) are optimized for each B decay mode, and depend on momentum. We form 
7]-0 _^ ryry caudidatcs from pairs of two photons. The photon energy in the laboratory frame 
is required to be greater than 50 MeV. The 7r° candidates should have invariant masses in 
the range 0.118 GeV/c^ < M^^ < 0.150 GeV/c^. 

We reconstruct a candidate by combining a —>■ ti^tt" and two tt^s (a Kg —>■ tt+tt^) 
for B^ — > KgTc^TT^ {KgKg) decay. Two kinematic variables are used for the B^ candidate 



selection: AE = Eb — -Ebeam and Mbc = y -f'bcam " Pb^ where -Ebeam is the cms beam energy, 
and Eb and pb are the cms energy and momentum, respectively, of a -B candidate. We 
require \AE\ < 0.25 (0.20) GeV and Mbc > 5.20 GeV/c^ for B^ K^-k^-k^ {K^s^^s)- 

We find 34% (< 1%) of events have multiple 5° ^ /^gTrOvrO (iT^/sT^) candidates. We 
choose the candidate having the smallest value for the 7r° for B^ K^t^^ti^, and having 
the smallest difference of M^+t^- from the nominal mass for B^ K^K^g. For B^ — >■ 
candidates, we find 23% of signal events in the Monte-Carlo (MC) simulation are 
incorrectly reconstructed self-cross-feed (SCF) events; for events at least one particle is 
replaced with one from either the accompanying B meson decay or from beam background. 
For 99.5% of the SCF events, only photons are replaced; the B vertex position is always 
reconstructed with the correct K^. Therefore, we do not distinguish between the SCF and 
correctly-reconstructed B candidates for this analysis. No SCF component is found for 
5° ^ KO/sfO decays. 

We remove B^ —>■ Kgn^rc^ candidates satisfying one of the following conditions: M^o^o > 
4.8 GeV/c^ to suppress 5° iT^vrO decays, 1.77 GeV/c^ < M^o^o < 1.94 GeV /J (D° 
veto), 3.27 GeV/c^ < M^o^o < 3.49 GeV/c^ (xco veto) and cos > 0.9, where 6^ is the 
angle between the 7r° boost direction from the laboratory frame and one of the photons in 
the 7r° rest frame. We impose the last condition to suppress 6 ^ 57 backgrounds. 

The dominant background is the continuum, e~^e~ qq {q = u, d, s, c). To suppress this 
background, we form a likelihood ratio 71 = Lsig/(Lsig + I^bg), where i^sig(bg) is a likelihood 
function for the signal (continuum) based on the event topology and the 5° flight direction 
with respect to the beam axis in the cms. The B'^ — > K'^tt^tt^ [K'^K'^) candidates having 
71 > 0.1(0.25) are selected. We also use 7Z to discriminate the signal and background in the 
fits. 

We identify the flavor of /tag using an algorithm |19|] that provides two variables: q defined 
in Eq. ([1]) and r. The parameter r ranges from r = for no flavor discrimination to r = 1 
for an unambiguous flavor determination. Candidate events having r > 0.1 are divided into 
six r-bins {i = 1,6). The wrong tag fraction in each i bin, Wi, and their differences between 
B^ and 5°, Awe, are determined using data [2o|. 



We apply the vertex reconstruction algorithm of Ref. [16|. The vertex position of the 
fcp decay is determined using the Kg — > tt+tt^ pseudo-track and the IP proflle. The two 
charged tracks are required to have enough associated SVD hits for vertex reconstruction. 
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We reconstruct the /tag decay vertex position using charged tracks that are not used 

for the fcp reconstruction. We measure both Sf and Af parameters with candidate events in 
which the vertex position is reconstructed successfully. Candidates with no vertex position 
information are used only for Af measurements. Using a signal MC sample, we estimate a 
vertex reconstruction efficiency of 45% (56%) for fC^TT^Tr^ (K^Kg) decays. 

We determine signal yields with an extended unbinned maximum likelihood fit that makes 
use of the variables AE, Mbc and 71. The selected candidate events include not only the 
signal but also continuum and BB backgrounds. 

We use a two-dimensional smoothed AE'-Mbc histogram for the B^ Kgn^Ti^ signal 
probability density function (PDF), while we model the AE (Mbc) shape using a sum of two 
Gaussians (a single Gaussian) for B^ — > KgKg decay. Histogram PDFs are employed for 
the TZ distributions of signal in both decay modes. The PDFs are determined using signal 
MC simulations. 

The AE (Mbc) shape of the continuum background is modeled as a polynomial (an 
ARGUS [2l|) function. The parameters of the background functions are fioated in the fit. 
An TZ PDF for continuum background is determined from a data sample collected at a center 
of mass energy 60 MeV below the T{4S) resonance (off- resonance) for the B^ K^vr^vr^ 
and events in the sideband region (Mbc < 5.26 GeV/c^) for B^ ^s-^s decay. We use 
the off-resonance data for the former decay mode, since BB background contaminates the 
sideband region, while the latter decay mode has a negligibly small contribution from BB 
background. A large MC sample of BB background is used to obtain a smoothed histogram 
AE-Mbc PDF and a histogram PDF of U for 5° K^tc^tc^ decay. 

By fitting the data, we determine the signal yields in the signal box. The signal box is 
defined as -0.15 GeV< AE < 0.1 GeV and Mbc > 5.27 GeV/c^ {\AE\ < 0.1 GeV and 
Mbc > 5.27 GeV/c^) for 5° ^ K^gir^n^ (^5^5) decay. For 5° ^ K^7r°7r° decay, we divide 
the candidate events into two categories: events with and without vertex information, and 
fit the data subsets separately to take into account a possible difference of the signal fraction. 
For B^ ^s-^s decay, instead of separating the candidate events in the fit, we perform 
a fit to the whole data sample and employ the vertex reconstruction efficiency estimated 
from a MC simulation (sideband events) for the signal (continuum) events to determine 
the signal fractions for the two categories. We obtain 129 ± 21 (178 ± 24) B^ KgH^ir^ 
events with (without) the vertex information, which is consistent with the expected vertex 
reconstruction efficiency, and 58 ±11 B^ ^ K^K^ events. Figure [T] shows the projections 
of AE, Mbc and TZ for candidate events. 

We determine the CP violation parameters from an unbinned maximum likelihood fit to 
the At distribution of the candidate events in the signal box. For each event, we calculate 
a likelihood value: 

Pi = foi) E fki^E, M^c, n)VkiAti) ® RkiAti) 

k 

+UVoi{AU), (2) 

where k indicates signal, continuum (qq) and BB background components, fk{AE, Mi,c,TZ) 
is the fraction of the component /c as a function of AE, M\yc,TZ. For B^ — * KgK^, the 
BB contribution is neglected, i.e. f^-^ = 0. For the signal function Vsig{At), we use 
Eq. ([1]) modified to incorporate the fiavor mis-assignment effect using we and Aw£. This 
distribution is convolved with a resolution function i?sig(At) [22] to obtain the signal PDF. 
The continuum At distribution Vqq{At) consists of prompt and lifetime components, which is 
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FIG. 1: (a) distributions of TZ for the AE and M^c signal region, (b) Mbc for the AE signal 
region with TZ > 0.9, and (c) AE for the Mbc signal region with TZ > 0.9 for — > Kgir^n^. 
(d) distributions of TZ for the AE and Mbc signal region, (e) Mbc for the AE signal region with 
TZ > 0.6, (f) AE for the Mbc signal region with TZ > 0.6 for B° K^K^- The sohd curves 
and histograms show the fits to signal plus background distributions, and hatched areas show the 
background contributions. 



convolved with a resolution function Rqg{At) comprised of two Gaussians; all the parameters 
of the functions are determined using the off-resonance (sideband) data for K^tt^tt^ 
{KgK'^) decay. We employ a large MC sample to obtain V^-^{At) and R^-^{At); the same 
functional forms as those of the continuum background are used. All the PDFs are combined 
with a outlier PDF Voi{At) that takes into account a small fraction (/oi) of events having 
large At values, giving a likelihood value Pj for i-th candidate event. For candidate events 
with no vertex information, we use the likelihood function given in Eq. ([2]) integrated over 
At. 

In the fit, Sf and Af are the only free parameters and are determined by maximizing 
£ = UiPi- The fit yields S^o^o^o = +0.43 ± 0.49, Ako^o^o = -0.17 ± 0.24, SkOk^ = 
—0.38 ± 0.77 and A^oj^o = —0.38 ± 0.38, where the errors are statistical. Figure [2] show the 
At distributions and the asymmetry Acp in each At bin, where Acp = (A^++A^_) 
and is the number of candidate events with q = +1 (—1). 

Table [I shows the systematic errors for the measured CP violation parameters. For the 
uncertainties on the wrong tag fractions, we vary the wrong tag fraction parameters, we and 
Aw^, by ±1(7 individually, and sum up the CP violation parameter variations in quadrature. 

The parameters r^o and Am^ are varied by their errors 23| to determine the systematic 
uncertainties from the physics parameters. 

We perform fits by varying each parameter in the signal resolution function i?sig(At). The 
differences between the fitted parameters and the nominal values are added in quadrature. 

We vary the background At PDF parameter values by ±lcr, and take the quadratic sum 
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-7.5 -5 -2.5 2.5 
At(ps) 



FIG. 2: At distributions and asymmetries for events with good tags (r > 0.5) for (a),(c) — > 
iCg7r°7r° with 7^ > 0.9 and (b), (d) B° K'^K^ with TZ > 0.6. In the At plots, the dashed 
and soUd Unes show the fit result for events with q = ±1. The dotted lines show the background 
contribution. In the asymmetry plots, solid curves show the fit results; the dashed curve in (c) 
shows the SM expectation ^ = and S = — sin2(/>i. 



TABLE I: Systematic errors of the measured CP violation parameters. 



sources 








•^R°R° 


wrong tag fraction 


±0.01 


±0.01 


±0.02 


±0.01 


physics parameters 


< 0.01 


< 0.01 


±0.01 


±0.01 


resolution function 


±0.04 


±0.02 


±0.06 


< 0.01 


background At shape 


±0.05 


±0.02 


±0.04 


±0.02 


background fraction 


±0.05 


±0.03 


±0.04 


±0.02 


possible fit bias 


±0.02 


±0.01 


±0.02 


±0.01 


vertex reconstruction 


±0.01 


±0.02 


±0.01 


±0.02 


tag side interference 


< 0.01 


±0.04 


< 0.01 


±0.03 


total 


±0.09 


±0.06 


±0.08 


±0.05 
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of the CP violation parameter variations as the systematic error. For — > i^^^vr'^Tr^, we 
also take into account the uncertainties originating from possible CP violation in the BB 
background. Using a large MC sample, we find 18% of the charmless B decay background are 
CP eigenstates. Each CP asymmetry parameter is varied between ±1, and the differences 
are added in quadrature for the systematic error. 

To estimate the background fraction systematic error, we vary each parameter in the 
PDFs of AE, Mbc and TZ, and add the differences in quadrature. For B^ — >■ i^gTr^vr", 
we also take into account uncertainties in the BB background yields; we vary the charm 
(charmless) B decay background yields by ±2o" (±100%). For B^ —>■ KgTc^Tc^, systematic 
uncertainties coming from the possible signal shape difference due to the uncertainty in the 
intermediate resonant state are also included. We generate pseudo-experiment MC samples 
based on PDFs assuming 100% resonant decays for the cases of B^ — > foK^ and B^ — > K*^7r^ 
decays, fit the samples with the nominal PDFs, and take the differences as systematic errors. 

A possible fit bias is examined by fitting a large number of MC events. We find no 
statistically significant bias and assign the MC statistical error as the systematic error. 

The systematic uncertainties for the vertex reconstruction are estimated by changing the 
charged track selection criteria, criteria on the vertex fit requirements, small bias correc- 
tions on Az, smearing due to the B flight length, and the SVD mis-alignment parameters. 

Finally, we take into account the possible CP violation effect in B^ /tag decay j23 |. 

To validate our measurements, various checks are carried out. We measure the B^ lifetime; 
we obtain tbo = 1.32 ± 0.27 (1.58 ± 0.44) ps for 5° K^Ti^n^ (^5^5) decay, in agreement 
with the world average (WA) value [23]. We perform a fit to sideband events; no asymmetries 
are found for both decay modes. We measure the branching fractions of B^ D^{^ 
ir°7r°)7r° and 5° ^ i^s^^^; we obtain (3.3 ±0.4) x 10"^ and (1.1 ±0.2) x 10"^ respectively, 
consistent with the WA values j^, [2^. The errors are statistical only. 

In summary we measure the CP violation parameters in FCNC 5° decays into K^tt^tt^ 
and KgK^ using a data sample containing 657 x 10^ BB pairs. We obtain 

Sko^o^o = ±0.43 ±0.49 ±0.09, 

Ako^o^o = -0.17 ±0.24 ±0.06, 

Skoko = -0.38 ±0.77 ±0.08 and, 

^^0^0 = -0.38 ±0.38 ±0.05, 

where the first and second errors are statistical and systematic, respectively. We estimate 
that the deviation from the SM expectation for CP violation parameters in 5° — i^^^vr^vr^ 
decay ^ = and S = sin2(/)i has a significance of 2. Oct. No CP asymmetry is found in 
B^ —>■ KgKg decay, which is consistent with the SM prediction. Our results are consistent 
with other measurements (7l.l9l.i0l. 
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